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ABSTRACT 



We present HST WFPC2 observations in three bands (F450W=B, F467M and 
F814W=I) of a group of three galaxies at z = 2.8 discovered in a ground-based narrow- 
ed \ band search for Lya emission near the z — 2.8 quasar PKS0528 — 250. One of the 

galaxies is a damped Lya (DLA) absorber and these observations bear on the relation 
between the DLA clouds and the Lyman-break galaxies and the stage in the evolution 
of galaxies they represent. We describe a procedure for combining the undersampled 
WFPC2 images pointed on a sub-pixel grid, which largely recovers the full sampling of 
the WFPC2 point spread function (psf). These three galaxies have similar properties 
to the Lyman-break galaxies except that they have strong Lya emission. The three 
■ galaxies are detected in all three bands, with average tub ~ 26, mj ~ 25. Two of 

qq I the galaxies are compact with intrinsic (i.e. after correcting for the effect of the psf) 

Q\ . half-light radii of ~ 0.1 arcsec (0.4/i _1 kpc, q Q — 0.5). The third galaxy comprises two 

similarly compact components separated by 0.3 arcsec. The HST images and a new 
ground-based Lya image of the field provide evidence that the three galaxies are more 
extended in the light of Lya than in the continuum. Combined with the evidence from 
the Lya line widths, previously measured, this suggests that we are measuring the size 
of the surface of last scattering of the escaping resonantly-scattered Lya photons. The 
measured impact parameters for this DLA galaxy (1.17 arcsec), for a second confirmed 
system, and for several candidates, provide a preliminary estimate of the cross-section- 
weighted mean radius of the DLA gas clouds at z ~ 3 of < 13/i _1 kpc, for q = 0.5. 
The true value is likely substantially smaller than this limit as DLA clouds at small 
r> | impact parameter are harder to detect. Given the observed sky covering factor of the 

^ • absorbers this implies that for q Q = 0.5 the space density of DLA clouds at these 

redshifts is more than five times the space density of spiral galaxies locally, with the 
actual ratio probably considerably greater. For q = 0.0 there is no evidence as yet 
that DLA clouds are more common than spiral galaxies locally. We summarise evidence 
that filamentary structures occur in the distribution of galaxies at high redshift. 

Key words: galaxies: formation - quasars: absorption lines - quasars: individual: 
PKS 0528 - 250 



1 INTRODUCTION 

By studying typical galaxies at high redshift, here z > 2, 
we can record the origins of no rmal galaxies such as our 
own Milky Way. Schmidt (1965) was the first to observe 
a high— redshift galaxy when he obtained a spectrum of 



* Based on observations with the NASA/ESA Hubble Space Tele- 
scope, and on observations collected at the European Southern 
Observatory, La Silla, Chile 



3C9, a radio— loud quasar of redshift z = 2.01. Normal 
star— forming high— redshift galaxies are about 1000 times 
fainter, and remained undiscovered until the present decade. 
A small number of candidates (some of which may be ac- 
tive galaxies) have been identified with 4— metre telescopes, 
through the detection of Lya emission (e.g. Steidel, Sargent, 
and Dickinson 1991, Lowenthal et al. 1991, M0ller and War- 
ren 1993, Pascarelle et al. 1996, Francis et al. 1995). Another 
approach has been to employ deep broad— band imaging to 
identify candidates by the expected Lyman— limit discon- 
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timiity in their spectra (Steidel and Hamilton 1992). The 
brightest Lyman— break galaxies have ran ~ 24 and spec- 
troscopic confirmation has only become feasible with the 
completion of the Keck 10— metre telescope. The report by 
Steidel et al. (1996) on the results of spectroscopic observa- 
tions of Lyman— break candidates really marked the begin- 
ning of the statistical study of high-redshift galaxies. They 
were able to confirm redshifts for 15 star— forming galaxies in 
the range 3.0 < z < 3.5. They discovered that Lyman— break 
galaxies generally show weak Lya emission, but that the 
rest— frame ultra— violet spectra may be recognised by the 
presence of absorption lines characteristic of the spectra of 
nearby star— forming galaxies. Madau et al. (1996) applied 
the same methodology to the Hubble Deep Field data to 
extend the results to fainter magnitudes, and their analy- 
sis currently provides the best lower limits to the integrated 
star formation rate in galaxies at z > 2. 

Complementary to the searches for starlight from 
high— redshift galaxies have been the analyses of absorption 
lines in the spectra of high— redshift quasars. These studies 
have provided measurements of the mass density of neutral 
hydrogen in the universe (e.g. Wolfe et al. 1986, Lanzetta et 
al. 1991), and the metallicity of the gas (e.g. Pettini et al 
1994), and how these quantities have changed with redshift. 
The analysis of Pei and Fall (1995) of the absorption-line 
data reconstructs the global history of star formation, gas 
consumption, and chemical enrichment, accounting in a self- 
consistent way for the effects of the progressive extinction of 
the background quasars due to dust as star formation pro- 
ceeds. The advantage of the absorption-line approach to the 
history of star formation is that it is global in nature as all 
the neutral gas at any redshift is directly observed. With the 
deep imaging studies it is necessary to extrapolate below the 
survey flux limit to measure the total star formation rate. 

The weakness of the global approach however is that 
it tells us nothing about the sizes and morphologies of the 
galaxies in which the star formation occurs. Deep imaging of 
DLA absorbers provides the missing information. In other 
words if we imaged the absorbers we could combine the data 
of faint galaxy surveys (luminosities, sizes, shapes) with the 
information on the gas obtained from the spectroscopic stud- 
ies (column densities, chemistry), allowing a more detailed 
comparison with theories of galaxy formation. 

Unlike for Mgll absorbers (e.g. Bergeron and Boisse 
1991, Steidel, Dickinson, and Persson 1994a) little progress 
has been made in such a programme of imaging of DLA 
systems. There are only two published unambiguous detec- 
tions, discussed below, as well as a small number of candi- 
date counterparts of DLA absorbers (at both low redshift — 
Steidel et al. 1994b, Le Brun et al. 1997 - and high-redshift 
— Aragon-Salamanca et al. 1996). In this paper we report 
on the results of 30 orbits of imaging observations with the 
Hubble Space Telescope (HST) of a group of three galax- 
ies at z = 2.81, named SI, S2, S3, one of which (SI) is a 
damped Lya absorber. The three galaxies lie in the field of 
the quasar PKS 0528 — 250 , and were detected in the light 
of Lya emission by M0ller and Warren (1993, hereafter Pa- 
per I). A detailed discussion of their nature was presented 

by Warren and M0ller (1996, hereafter Paper II Q). The 
t We note here that recent spectroscopic observations by Ge et 



Table 1. Journal of HST observations 



Date 


Filter 


Band 


Exposure 








(sec) 


1995 Feb 2 


F467M 


M 


12900 


1995 Mar 5 


F467M 


M 


12900 


1995 Mar 8 


F814W 


I 


12900 


1995 Mar 8 


F450W 


B 


12900 


1995 Mar 9 


F467M 


M 


12900 



HST images provide measures of their sizes, magnitudes, 
and colours. The other high-redshift damped system that 
has been successfully imaged is the absorber at z — 3.150, 
t owar ds the quasar 2231 + 131, observed by Djorgovski et al. 
( 1996) ) . In this paper we intercompare the measured prop- 
erties of these two absorbers, the two companion galaxies 
S2 and S3, and the population of Lyman— break galaxies, to 
draw conclusions about the space density of DLA absorbers, 
their structure, and the relation between the DLA absorbers 
and the Lyman-break galaxies. 

The layout of the rest of the paper is as follows. 
The HST observations of the field towards PKS 0528 - 
250 are described in Section 2, as well as new ground— based 
narrow— band observations of the field. The three galaxies 
are very small, and the HST images were dithered with a 
half-integer pixel step in order to improve the image sam- 
pling. The algorithm for combining the images is outlined 
in this section. In Section 3 we present the reduced HST 
images. In Section 4 we provide the results of aperture pho- 
tometry and profile fitting of the HST and ground-based 
images. Finally, in Section 5 we discuss these results and 
their implications for our understanding of the nature of 
DLA absorbers and Lyman-break galaxies. 



2 OBSERVATIONS 
2.1 HST observations 

The HST observations were made in 1995 February and 
March, using the Wide Field Camera of WFPC2. The jour- 
nal of observations appears in Table 1. We used three differ- 
ent filters: six orbits each of observations with the WFPC2 
standard broad B (F450W) and I (F814W) filters, and 18 
orbits with the F467M medium-passband filter. The last fil- 
ter, hereafter M, has a width of 215A (FWHM), and is the 
narrowest WFPC2 filter that contains the wavelength of red- 
shifted Lya of the targets, 4630A. 

The 30 orbits were scheduled as five visits, each of six 
orbits. For each visit we observed through only one filter. 
In each case we obtained a 1900 sec exposure in the first 
orbit, and an uninterrupted 2200 sec exposure in each of 
the following 5 orbits. The six exposures of each visit were 
grouped in three pairs. The two exposures of each pair were 
obtained with exactly the same pointing. Relative to the 

al. (1997) confirm our conclusion that the Lya emission from the 
DLA absorber is due to star formation rather than photoionisa- 
tion by the quasar. 
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first, the second and third pairs were shifted, respectively, 
by 0.45 arcsec in x, and by 0.45 arcsec in both x and y. 
The step of 0.45 arcsec is exactly 4.5 WFC pixels, and this 
strategy was chosen with a view to improving the rejection 
of cosmic rays, and to allow an increase in the sampling, 
since the WFPC2 point spread function is undersampled by 
the 0.1 arcsec WFC pixels. The latter point is important be- 
cause we need good sampling for a morphological study, but 
our objects are too faint to be observed with the Planetary 
Camera. 

For the basic reduction of the data we obtained, from 
the STScI data archive, the best biases, darks and flats ap- 
propriate for the observing dates and applied these in stan- 
dard ways. We also obtained the relevant hot pixel lists, 
and created a bad pixel mask and a hot pixel mask for each 
image. 

To combine the frames for each filter we used the 
weighting scheme and the sigma clipping algorithm de- 
scribed in Paper I, which is optimal for faint sources. We 
have extended the algorithm to allow for half integer shifts 
between input frames. A short description of the procedure 
follows. A detailed discussion will be provided elsewhere 
(M0ller, in preparation). 

The sub— pixel combination is undertaken in two stages, 
of which the second is iterative. The first stage is the com- 
bination of the frames using integer pixel shifts, with the 
rejection of cosmic rays by sigma clipping, and the mask- 
ing of bad pixels. In this combined image each pixel is then 
divided two by two into sub— pixels, providing the first es- 
timate of the final sub— pixelised image. The pixels in the 
original data frames are then similarly divided, and each 
frame is registered with the combined image, to the nearest 
sub— pixel. The sub-pixel registration allows improved re- 
jection of cosmic rays in the subsequent processing, notably 
where there are strong gradients in the data e.g. in images 
of stars. 

The first combined image provides a prescription for 
each original data pixel as to the distribution of the de- 
tected incident photons within the pixel. In each frame, in 
each pixel, the counts are then distributed to the four corre- 
sponding sub— pixels according to this prescription. A sec- 
ond combined frame is then constructed, again after appro- 
priate rejection of bad pixels and cosmic rays. This new 
frame provides an improved prescription for the redistri- 
bution of counts within pixels, and the final frame is con- 
structed in this manner through iteration. The iteration is 
halted as soon as the mean x 2 per pixel has reached a sta- 
ble minimum, typically after six to twelve iterations. Along 
with the combined image, the corresponding variance image 
is created. 

The measured profiles of objects in the final frame will 
be similar to the profiles that would have been obtained for 
the same instrument but with pixels of side half the size. 
Thus the algorithm improves the sampling by a factor of 
two, but does not attempt to deconvolve the resolution set 
by the optics. It does, however, remove the extra "box-car" 
smoothing of the image caused by the finite pixel size. The 
algorithm is strictly locally and globally count conserving, 
and the observed number of counts recorded in a whole de- 
tector pixel is never changed (as opposed to deconvolution 
algorithms attempting to correct for the effect of the psf). 

Most cameras used in optical astronomy have pixels 
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Table 2. Journal of NTT/SUSI observations 



Date 


Filter 


Exposure 






(sec) 


1993 Dec 10/11 


N 


15700 


1994 Nov 30/Dec 1 


N 


18000 


1994 Dec 1/Dec 2 


N 


18500 



which are small compared to the resolution FWHM of the 
images, and hence their intrinsic pixel smoothing is insignif- 
icant. For WFPC2 the pixels are of the same size as the 
resolution, so the pixel smoothing adds significantly to the 
final resolution if it is not corrected for. The correction al- 
gorithm outlined above therefore, for WFPC2, has the ad- 
vantage of improving both sampling and resolution, but it 
does not suffer the disadvantage of mixing information be- 
tween unrelated detector pixels as does psf-deconvolution 
algorithms. The main difference between our algorithm and 
that used by the HDF group (Hook and Fruchter, 1997) is 
that the HDF group decided not to correct for the pixel 
smoothing because the correction introduces a small corre- 
lation of noise on the scale of an original pixel. The HDF 
images therefore have somewhat lower resolution but statis- 
tically independent sub-pixels. 

2.2 NTT observations 

We obtained a new deep ground-based Lya image of the 
field using the SUSI instrument on the ESO New Technol- 
ogy Telescope (NTT) , over two runs in December 1993 and 
December 1994. The journal of the observations appears in 
Table 2. The narrow-band filter used for the NTT observa- 
tions is the same filter, hereafter N, as was used for the ESO 
3.6m observations reported in Paper I, and has a central 
wavelength A = 4633 A, FWHM 23 A. 

Our original ESO 3.6m narrow-band image was badly 
affected by filter ghosts, and by the poor seeing, average 1.7 
arcsec FWHM. The purpose of the SUSI image was to ob- 
tain a deeper ghost-free image (as the filter is not inclined 
in SUSI) with better seeing. The observations were broken 
into 1 hour exposures obtained with small shifts in point- 
ing. The data reduction and the combination of frames fol- 
lowed the procedures described in Paper I. The images in 
the final narrow-band frame have FWHM 0.96 arcsec, with 
a sampling of 0.26 arcsec per pixel. The source SI lies at 
an angular separation of only ~ 1.2 arcsec from the line of 
sight to the quasar. Due to the relatively good seeing the 
faint image of the quasar (from the small amount of quasar 
light that leaks through the filter wings) and the image of SI 
are separate. There are several unsaturated images of bright 
stars in the SUSI frames, so that using DAOPHOT we were 
able to subtract the image of the quasar. 

3 HST IMAGES 

In Figure 1 we show a small section of the field around the 
quasar, that includes the three galaxies SI, S2, S3. The fig- 
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4.1 Aperture photometry 



See attached GIF colour image. 



Figure 2. HST Wide Field Camera sub-pixelised images showing 
the quasar and the damped Lya galaxy SI. Each box is 3.5 arcsec 
on a side, and the size of each sub-pixel is 0.05 arcsec. Top row 
combined images for, from left to right, the B, M, and I niters. 
Bottom row, the same images after subtraction of the quasar psf, 
to disclose the image of SI. 



ure is a weighted sum of the final HST frames for the three 
filters. The three galaxies are clearly visible in this figure, 
and are also each detected in the individual summed frames 
for each filter. The galaxy SI lies at an angular separation 
of only 1.17 ± 0.02 arcsec from the quasar, and is only vis- 
ible after subtraction of the quasar image. In Figure 1 we 
have inserted a small panel from the image with the quasar 
subtracted, to reveal SI. The lower panel of Figure 1 is the 
same image lightly smoothed to enhance the contrast of faint 
objects. 

For the psf subtraction we used Tiny TIM v4.1 model 
psfs. Model psfs were made for each filter, and were con- 
structed on the 0.05 arcsec sub-pixel scale. Images of SI in 
all three bands, before and after psf subtraction, are shown 
in Figure 2. Construction of sub-pixelised images and psfs 
provides a substantial improvement in the psf subtraction 
over our preliminary integral pixel shift reductions shown in 
M0ller and Warren (1996). 

A comparison between images obtained in different 
bands is provided in Figure 3, which shows close-up images 
of all three galaxies. The sources SI and S2 are both com- 
pact, without any obvious substructure, while the source S3 
is composed principally of two compact components, sep- 
arated by 0.3 arcsec. A possible third component is just 
visible, but contains only 4% of the light, and has not been 
included in the profile fitting described in the next section. 

At this point it is worth commenting on the low signal- 
to- noise ratio (S/N) apparent in Figure 3 of the M-band 
images of sources S2 and S3 compared with the B-band and 
I-band images. The Lya emission line acounts for a large 
fraction of the M-band flux for the three sources. As detailed 
in Section 5 we find evidence that the Lya emission is more 
extended than the continuum emission for the three sources, 
which would explain the low S/N of the M-band images. 



4-1.1 Aperture photometry of HST images 

We have measured magnitudes for the three galaxies us- 
ing both a small (isophotal) aperture, to obtain accurate 
colours, and a large (circular) aperture, to measure total 
fluxes. The results are provided in Table 3, as well as the 
colours tub — mi, tub — rriM , for each source, for both aper- 
tures. The magnitudes are on the HST system, with zero 
points taken from Holtzman et al. (1995) (their Table 9). 

For the small— aperture measurements we selected 
isophotes that defined areas of 0.13, 0.13, and 0.31 square 
arcsec for the galaxies SI, S2, S3, respectively. For each ob- 
ject the identical aperture was applied in each passband. The 
quoted photometric errors account for photon noise, read- 
out-noise, and the error associated with the uncertainty in 
the determination of the local sky level. 

The large aperture used was a circle of diameter one 
arcsec. The zero points of the HST photometric system are 
defined such that for point sources aperture magnitudes in 
a circular aperture of this size provide total magnitudes. 
Therefore our large-aperture magnitudes will be total mag- 
nitudes if the sources are compact. This is certainly the case 
for the continuum emission, as demonstrated by the profile 
fits described in the next subsection. The measurement of 
the flux within the large aperture is difficult for SI. This 
is because the aperture includes regions where the subtrac- 
tion of the quasar profile is very uncertain. Instead, for this 
source we applied the aperture to the galaxy-profile models 
described in the next subsection, but have not attempted to 
quantify the photometric errors. In computing the best— fit 
models the bad regions were ignored. Therefore, for SI, the 
models ought to provide a more reliable estimate of the flux 
within the larger aperture. 



4-1-2 Aperture photometry of NTT image 

For the NTT image we used a circular aperture of diam- 
eter 2.6 arcsec. We measured the aperture correction be- 
tween this diameter and a very large diameter for bright 
stars in the frame, and applied the correction to the mea- 
sured magnitudes of our sources. Again, then, if the sources 
are compact this procedure provides total magnitudes. The 
photometry for the NTT observations is on the AB system. 

Also provided in Table 3 are the values of the rest-frame 
equivalent width (EW) of the Lya emission line for the three 
sources. These were computed by creating power-law spectra 
/„ oc is , with emission lines of different EW, and applying 
absorption at wavelengths shortward of the Lya line centre 
as approporiate for z = 2.81, using the models of M0ller 
and Warren (1991). The line EW was adjusted to reproduce 
the measured large-aperture ms — mjv colours. The value 
quoted in the table is the model unabsorbed EW. We used 
the same model spectra to compute the expected value of 
the colour ms — mM, quoted in the last row of Table 3. These 
predicted colours are used in Section 5.2, in considering the 
evidence that the Lya emission is more extended than the 
continuum emission. 



4 PHOTOMETRY 



4.2 Galaxy profile fits 
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See attached GIF colour image. 



Figure 1. HST Wide Field Camera 2 sub— pixelised image of a region 43arcsec X 17arcsec around the quasar PKS0528 — 250 (marked 
Q), showing the three galaxies SI, S2, S3. The size of each sub-pixel is 0.05 arcsec. The upper panel is a weighted average of the data for 
the three filters, from all 30 orbits. The lower panel is the same image convolved with a Gaussian profile of a = 1.3 sub-pixel. An inset 
panel shows a small region to the S of the quasar after subtraction of the quasar image, revealing the image of SI. The scale is shown by 
a horizontal bar of length 3 arcsec, the orientation is given by the N/E arrows to the right. 



Table 3. HST and NTT photometry of SI, S2 and S3 

Small (isophotal) aperture Large (circular) aperture 

Filter SI S2 S3 Sl a S2 S3 



HST F450W 


m B 


25.90+q^ 


26 70+ 09 
( U — 0.08 


25.63l°;° 5 5 


25.45 


26.281°;^ 


25.46l°;°« 


HST F467M 


m M 


25.56±S;S 




25.43lg;f 2 


25.18 


25.76l°-g 


25.161°;" 


HST F814W 


mi 


24.96+° 


25.961°;°? 


24 94+O.O6 


24.85 


25.66l°« 


24.881°;°? 




m B - 


-mi 


94+°' 30 


74+ - 11 


0.69l°;g 7 r 


0.60 


62+ ' 20 


0.581°;!! 




m B ~ 


-m M 


S4+°' 27 


0.08t°|2 


0.20l°;« 


0.27 


0-52±g:i| 


0-30l°i? 


NTT N 


m 


N 








23.071°;!° 


23.481°;" 


23 17+ - 11 




m B - 


-m N 








2.38 


2.801°-, 


2 29+ ' 14 


EW rost 












51 


88i 2 21 


46l° 


Expected 


rn B - 


-m M 








0.56 


0.77l°; 1 1 2 


52+ - 06 



a B, M, I values measured from the model galaxy profiles 
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See attached GIF colour image. 



Figure 3. This figure presents a comparison of images of SI (top row), S2 (middle row) and S3 (bottom row) in all HST bands. From 
left to right the first three columns show the B, M, and I filter images. For the source SI the quasar has been subtracted, and the low 
S/N F467M images have been lightly smoothed. In column 4 we show a weighted sum of the images from all three filters, in column 5 
our best model fit to the I band image, and in column 6 the residuals from this fit. All images in this figure are 2.0 by 2.0 arcsec with 
orientation as in Figure 1. 



Table 4. Galaxy profile fits 



object 


band 


r .5 


^0.5 


ro.s 


r .5 






dc Vauc. 


exp. 


average 


average 






arcsec 


arcsec 


arcsec 


kpc 


SI 


B 


0.21 


0.10 






SI 


M 


0.20 


0.09 


0.13 6 


0.48ft- 1 


SI 


I 


0.13 


0.08 


±0.06 


±0.21ft- 1 


S2 


B 


0.10 


0.06 


0.08 


0.29ft- 1 


S2 


I 


0.08 


0.07 


±0.01 


±0.05ft- 1 


S3-a 


B 


0.07 


0.06 


0.06 


0.23ft- 1 


S3-a 


I 


0.07 


0.06 


±0.01 


±0.03ft- 1 


S3-b 


B 


0.13 


0.09 


0.11 


0.42ft- 1 


S3-b 


I 


0.13 


0.10 


±0.02 


±0.08ft- 1 



a h = Ho/100, q a = 0.5 

6 average of values for B and I bands only 



4-2.1 Profile fits for HST images 

We have measured the half-light radii r 5 of the galaxies SI 
(after subtraction of the quasar image) and S2, and of the 
two components of the galaxy S3 by fitting de Vaucouleurs 
and exponential profile functions. The results are provided in 
Table 4. The two components of S3 are separated by 0.30 ± 
0.01 arcsec. The S/N of the M-band images of S2 and of S3 



were too low to provide meaningful results. There follows a 
description of the fitting procedure. 

Firstly any residual sky counts (or large-scale psf resid- 
uals for SI) in the frames were subtracted by means of a 
large-scale running modal filter, where the mode is estimated 
by fitting a Gaussian to the peak in the histogram of the 
counts in the sliding box. A box of side 51 sub-pixels (2.55 
arcsec) was used for S2 and S3. A smaller box of side 41 sub- 
pixels was used for SI, ignoring pixels close to the quasar 
centroid. For the profile fits themselves boxes of side 21 sub- 
pixels were used for S2 and S3, and of side 15 sub-pixels for 
SI. The results were fairly insensitive to the choices of box 
size. The model profiles are characterised by six parameters: 
x, y, ellipticity, position angle, r .s, and surface brightness at 
r .5- A mask was applied to exclude certain pixels from the 
fit, as required; for example the faint third source in fitting 
S3, and in fitting SI the region over which the quasar im- 
age is bright (i.e. where the Poisson noise is high, and where 
large residuals from the psf subtraction remain) . Trial model 
galaxy light profiles were convolved with the normalised sub- 
pixelised light profile of a nearby star, and the minimum \ 2 
fit was found for the six free parameters. Examples of the 
best— fit galaxy models are shown in Figure 3, as well as the 
images after subtraction of the models. 

Referring to Table 4 it can be seen that for each galaxy, 
for a particular profile function, there is good agreement in 
the measured half-light radii for the different filters, indicat- 
ing that the point spread function is adequately sampled in 
the sub-pixelised images. The goodness of fit of both pro- 
file functions is in all cases similar, so there is no indication 
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that either the de Vaucouleurs or exponential function is pre- 
ferred. In all cases a larger value of r .s is measured for the 
de Vaucouleurs function, and the discrepancy is larger than 
the scatter for the same function for different filters. This is 
particularly noticeable for galaxy SI and it is possible that 
this is caused by residuals left from imperfect subtraction 
of the quasar image. The purpose of fitting de Vaucouleurs 
and exponential functions is to be able to correct for the 
psf. The profile fitting indicates that the dominant uncer- 
tainty in the measurement of the intrinsic half-light radius 
is the choice of profile rather than Poisson noise. Therefore 
the range of half-light radii quoted in Table 4 for the differ- 
ent profile functions provides a measure of how accurately 
the half-light radius is determined. For this reason we have 
chosen to average all the measurements as the best estimate 
of the half-light radius, and to quote the scatter amongst 
the measurements for each object as a measure of the un- 
certainty, as listed in column 5 of the table. 

The psf has a significant effect on the measured half- 
light radii for the smallest objects. For example the half-light 
radius for S2 measured from the images, without correction 
for the psf, is 0.15 arcsec, which is about double the intrinsic 
value. 



4.2.2 Profile fits for NTT image 

The N-band ground-based images are of higher S/N but 
lower spatial resolution than the HST M-band data. Be- 
cause the continuum contributes only a small fraction of the 
light to the N-band images the N-band data essentially pro- 
vide the Lya profiles of the three sources, and we attempted 
to measure the profiles following the fitting procedure de- 
scribed above. For S2 and S3 the results are not particularly 
useful as evidence of whether the Lya emission is more ex- 
tended than the continuum emission. For both sources the 
measured intrinsic sizes are consistent with the continuum 
sizes measured from the HST images, but also consistent 
with the hypothesis that they are more extended. The N- 
band image of SI on the other hand is clearly extended. We 
measured half-light radii of 0.5li^o8 arcsec and 0.661° '^g 
arcsec for the exponential and de Vaucouleurs profiles re- 
spectively. The quoted errors are computed as follows: the 
value of the radius is varied in steps away from the best-fit 
value, and for each fixed value of the radius a new mini- 
mum x 2 fit is computed by allowing the other parameters 
to vary. The radius at which the \ 2 is 10 greater than the 
best-fit x 2 is identified. The change in radius relative to the 
best-fit value is the la error. For each object the results of 
the profile fitting demonstrate that the 2.6 arcsec diameter 
aperture used for the photometry includes most of the light 
from the source, so we will treat the measurements as total 
magnitudes. 



5 DISCUSSION: DLA ABSORBERS AND 
LYMAN-BREAK GALAXIES 

5.1 Continuum emission 

The three galaxies SI, S2, S3 are similar in their properties 
to members of the recently detected population of Lyman- 
break galaxies, except that they have strong Lya emission, 



with mean restframe EW = 60A. For example the B and I 
magnitudes for our three sources thb ~ 26, mi ~ 25, are 
within the range of magnitudes of the Lyman-break galax- 
ies of redshift z ~ 3 observed by Steidel et al. (1996) and 
Lowenthal et al. (1997). The average half-light radius of the 
four objects listed in Table 4 is 0.1 arcsec. This is smaller 
than the value quoted for the Lyman-break galaxies (Gi- 
avalisco et al. 1996, Lowenthal et al. 1997), by a factor of 
two to three. However we have broken S3 into two sub com- 
ponents and measured the radii for each subcomponent. In 
addition the radii quoted in Table 4 have been corrected for 
the effect of the psf. 

Young star-forming galaxies should have approximately 
constant flux per unit frequency over the restframe ultravio- 
let and optical regions of the spectrum. The expected colour 
of a source at z = 2.8 with a power-law continuum vary- 
ing as f v oc v°, after accounting for typical absorption in 
the Lya forest, is tub — mi = 0.66. As seen in Table 3 the 
niB — mi colours for the three sources SI, S2, S3 are consis- 
tent with this value, for both the small and large apertures. 

The fact that SI has properties similar to the Lyman- 
break galaxies suggests that DLA absorbers and Lyman- 
break galaxies may commonly be associated with each other. 
Lowenthal et al (1997) have also made a connection be- 
tween the two populations by suggesting that the broad 
Lya absorption line seen in several cases in the spectra of 
Lyman-break galaxies is a damped line from neutral hydro- 
gen in the galaxy. However the absorption lines seen in their 
average spectrum and in the relatively high S/N spectrum 
presented by Ebbels et al. (1996) are not optically thick in 
the line centre. The absorption-line profiles are difficult to 
interpret because the background source is extended. It is 
even possible that in some cases the absorption line is stel- 
lar absorption in the spectra of late B stars and that the 
galaxies are observed in a post-burst phase, rather than ac- 
tively forming stars (see Valls-Gabaud 1993) i.e. for some 
cases this may be the explanation for the weakness of the 
Lya emission line in the Lyman-break galaxies, as opposed 
to dust. High-redshift galaxies with strong Lya emission, like 
the three sources near PKS 0528 — 250 , may therefore be 
younger than many of the Lyman-break galaxies. 



5.2 Extended Lya emission 

In this subsection we summarise the evidence that the 
Lya emission from the three sources is more extended than 
the continuum emission. 

The evidence for S2 and S3 comes from an examina- 
tion of the m_B — mm and ms — tun colours, summarised in 
Table 3. The bottom row in Table 3 provides the predicted 
ms — m,M colours for each source computed, as detailed 
in Section 4.1.2, from the B and N total magnitudes. The 
small aperture ms — mu colours are smaller than the pre- 
dicted values, as would be the case if the Lya emission is 
more extended than the continuum. For S2 the difference 
is significant at 2.8er, and for S3 at 2.3(7. Notice also that 
the large aperture ms — tom colours are larger than the 
small aperture colours, which is consistent with the hypoth- 
esis of extended Lya emission. A broad absorption line in 
the M passband could only explain part of the discrepancy 
between the measured and predicted colours. For example 
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and size differences quoted above are only marginally signif- 
icant, these results accord with our earlier conclusion (Pa- 
per II) that the cause of the relatively broad Lya emission 
lines for SI and S3 is resonant scattering, as the photons 
escape though a high column density of HI. In this pic- 
ture the Lya image records the last scattering surface, which 
will be larger than the region of star formation. We note 
in passing that this implies that the observed strength of 
Lya absorption and emission in spectra can depend on the 
slit width used. 



See attached GIF colour image. 



Figure 4. Illustration of the evidence for extended Lya emission 
for, from top to bottom, SI, S2, S3. Each panel shows the same 
2 arcsec by 2 arcsec region as Figure 3, here smoothed by a 0.15 
by 0.15 arcsec boxcar filter to enhance faint features. The left 
column shows the summed image from all three filters. The middle 
column shows the M-band images with the computed contribution 
of the continuum subtracted. This image, then, contains only the 
Lya line component. In the right column we show how the Lya- 
only image would have appeared if the Lya profile were identical 
to the continuum profile. 

for S2, simply removing 50A of the continuum flux within 
the M band changes the tub — ttim colour by only 0.1 mag. 

The total B magnitude for SI (i.e. the large-aperture 
measurement from the HST image, Table 4) is not reliable 
because of the uncertainty associated with the psf subtrac- 
tion. This prevents us from usefully making the same com- 
parison between predicted and measured ms — tum colours 
for SI. However there is an indication that SI may also be 
more extended in Lya . The measurements of the half-light 
radius of SI from the ground-based Lya image yielded values 
of 0.5liQ gg arcsec for the exponential profile and 0.66io' 2 g 
arcsec for the de Vaucouleurs profile (Section 4.2.2). These 
values are significantly larger than the half light radius of 
the region of continuum emission of 0.13±0.06 arcsec (Table 
4), measured from the HST image. Such a diffuse contribu- 
tion of Lya to the M-band flux would be hard to detect, and 
could have been partially subtracted in attempting to re- 
move large-scale residuals from the psf subtraction (Section 
4.2.1). However we cannot rule out the possibility that there 
is also a similar diffuse contribution to the continuum flux. 

The argument for extended Lya is made visually in Fig- 
ure 4. Here we compare the M-band images after subtraction 
of the computed contribution of the continuum (i.e. show- 
ing the Lya emission only), against how the objects would 
have appeared if the Lya light profile were the same as the 
continuum profile. 

To summarise, there is evidence that the regions of 
Lya emission from these three sources are more extended 
than the regions of continuum emission. Although the colour 



5.3 The sizes of the DLA gas clouds 

In this section we make a summary of the measured impact 
parameters b of DLA absorbers. Using this we estimate the 
cross-section- weighted mean radius Rdla of the gas clouds 
at z > 2. Combining this information with the measured 
line density of DLA absorbers dn/dz we are able to infer 
the ratio of the comoving space density of DLA absorbers 
at z > 2 to the local space density of spiral galaxies. 

There have been many searches for optical counterparts 
of the DLA absorbers. In Table 6 we have listed the small 
number of DLA absorbers, with redshifts z > 0.8, for which 
likely candidate or confirmed optical counterparts have been 
detected. Listed are the quasar name, the absorber redshift 
and column density, the measured impact parameter (i.e. 
the projected physical separation between the optical coun- 
terpart and the line of sight to the quasar), and whether or 
not the counterpart has been confirmed, i.e. the redshift of 
the counterpart has been measured to be the same as the 
redshift of the absorber. 

In Figure 6 we plot impact parameter against column 
density for the absorbers listed in Table 6, for q = 0.5. The 
detection of the counterpart to a DLA absorber by broad- 
band imaging, which requires the digital subtraction of the 
quasar image, is more difficult for small impact parameters 
where the photon noise from the quasar image is greater. 
Therefore it is very likely that the measured impact pa- 
rameters of the few absorbers for which counterparts have 
been detected are larger than the average impact param- 
eter for DLA absorbers. The vertical line in Figure 6 at 
log(N(HI)) = 20.3 marks the conventional lower limit of the 
column density of DLA absorbers. The distribution of points 
in the figure is consistent with the expected anticorrelation 
between impact parameter and column density. We are in- 
terested in the mean impact parameter \)dla that would 
be measured if counterparts were detected for all DLA ab- 
sorbers in a large unbiased survey. From inspection of Fig- 
ure 6 we suggest that it is safe to conclude that for 2 > 2, 
q — 0.5, boLA is less than 7h~ kpc. For q = 0.0 we suggest 
boLA < 13/jT 1 kpc. The actual mean impact parameters are 
probably substantially smaller than these limits. 

We relate the mean impact parameter to the cross- 
section-weighted mean radius Rdla of the absorbers by 
bo la = olRdla- The edge of the DLA is the point at which 
the column density falls below N(HI) = 2 x 10 2() cm" 2 , 
and we effectively assume that the column density falls off 
sharply at larger radii. For face on disks a = |, and for a 
disk that is nearly edge on a = so we take a = 0.55 as 
an average value for randomly inclined disks. On this basis 
we infer the following limits on the cross-section-weighted 
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Table 5. Measured impact parameters b of DLA absorbers 
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Zabs 


log(N(HI)) 


b 


b 


confirmed 


reference 








{qo = 0.0) 


(qo = 0.5) 










enr 2 


kpc 


kpc 






0454 + 0393 


0.86 


20.8 


4.1ft 1 


3.3ft- 1 


N 


1 


0302 - 223 


1.01 


< 20.0 


6.2ft- 1 


4.9ft- 1 


N 


1 


1331 + 170 


1.78 


21.2 


4.7ft- 1 


3.1ft- 1 


N 


1 


0151 + 048A 


1.93 


20.4 


7.7ft- 1 


5.0ft- 1 


Y 


2 


1215 + 333 


2.00 


21.0 


8.4ft- 1 


5.3ft- 1 


N 


3 


0841 + 129 a 


2.37 


21.3 


8.0ft- 1 


4.7ft- 1 


N 


3,4 


0841 + 129 a 


2.48 


21.0 


8.0ft- 1 


4.7ft- 1 


N 


3,4 


0528 - 250 6 


2.81 


21.3 


8.1ft- 1 


4.5ft- 1 


Y 


5 


2231 + 131 


3.15 


20.0 


15.7ft- 1 


8.2ft- 1 


Y 


5 



° the detected galaxy could be the counterpart to one or other of the two DLA absorbers listed 
b impact parameter taken from HST image (this paper) 

References: 1. Le Brun ct al. (1997), 2. Fynbo ct al. (1997), 3. Aragon-Salamanca ct al. (1996), 4. Wolfe ct al. (1995), 5. M0llcr and 
Warren (1993), 6. Djorgovski ct al. (1996). 



mean radius of DLA absorbers at high redshift: 
Rdla < 23.6/T 1 kpc (z > 2, q a = 0.0) 

Rdla < 12.7/i _1 kpc (z > 2,q a = 0.5). 

Using these limits and the measured line density of ab- 
sorbers dn/dz we can compute limits to the ratio of the 
comoving space density of DLA absorbers to the the local 
space density of spiral galaxies. We follow essentially the 
methodology employed by Wolfe et al. (1986), and Lanzetta 
et al. (1991, hereafter LTW). For spiral galaxies locally they 
adopted a galaxy luminosity function $(7^) of Schechter 
form, with power-law index s, a power-law (Holmberg) re- 
lation between radius and luminosity, of index t, and a ra- 
tio between gas radius and optical (Holmberg) radius £, in- 
dependent of luminosity. They found that the incidence of 
DLA absorbers per unit redshift dn/dz at 2 ~ 2.5 was con- 
siderably higher than expected, by a factor F ~ 5, on the 
basis of no evolution in galaxy cross section or luminosity 
function normalisation. We now allow for evolution by sup- 
posing that the space density of DLA absorbers is higher 
than the local space density of spirals by a factor E<s>(z), 
and that the gas radii of galaxies are larger at high red- 
shift by a factor E r (z). Since dn/dz is proportional to the 
product of the space density and the galaxy cross section a, 
we have that F = E<s,(z)Er(z). By comparing the expected 
value of the cross-section-weighted radius R to the measured 
limits to Rdla, we obtain limits to E r (z). Then from the 
measured values of F we determine limits to E<s>(z), which 
is our goal. 

Under the above assumptions the cross-section- 
weighted average radius of DLA absorbers is given by: 

[z) Jo" 

_ E r (z)£,R,T(l+3t - s) 
~ T(l + 2t-s) 

where R* is the optical radius of a local i» spiral galaxy. 

Following LTW we adopt the following values of the pa- 
rameters: t = 0.4, s = 1.25, £ = 1.5, R* = 11.5/z -1 kpc. This 



leads to R(z) — ll.Oh -1 E r (z). Comparing with the above 
measured limits to Rdla we obtain the following limits to 
the growth factor of galaxy disks: 

E r (z) < 2.15 (2 > 2, go = 0.0) 

E r (z) < 1.16 (2 > 2, g = 0.5) 

For their sample D2, of which at least 30 out of the 38 
candidate DLA absorbers have been confirmed, LTW found 
F = 3.8, for g D = 0.0, and F = 7.1, for g D = 0.5. These re- 
sults then imply that the ratio of the comoving space density 
of DLA absorbers at high redshift to the local space density 
of spiral galaxies is given by: 

E*(z) = (D L A) /$, (spiral) > 0.8 (2 > 2,g D = 0.0) 

E*(z) = <f>*(DLA)/<f>*(spiral) > 5 (2 > 2,g D = 0.5) 

Because the actual average impact parameters of DLA ab- 
sorbers are very likely substantially smaller than the quoted 
limits, the actual comoving space density ratios are probably 
considerably greater than the limits quoted above. 

In Figure 5 we show also the relation between im- 
pact parameter and column density measured by Katz et 
al (1996) from a hydrodynamic simulation of a CDM uni- 
verse with g = 0.5. Given the limited spatial resolution of 
the simulation, and the fact that star formation and conse- 
quent feedback were not treated, it would be premature to 
draw any conclusions about the apparent good agreement 
between the results of the simulation and the observations. 
Nevertheless this plot and the above calculation underscore 
the importance of measuring impact parameters for a large 
sample of DLA absorbers, a goal we are pursuing by imaging 
with the STIS instrument on HST. 

5.4 The structure of the DLA absorbers 

The observations summarised above provide some indica- 
tion of the typical structure of a DLA absorber. The fol- 
lowing sketch is suggested. At the centre of the gas cloud, 
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Figure 5. Relation between impact parameter and column den- 
sity for all confirmed or candidate counterparts of DLA absorbers 
of redshift z > 0.8 (listed in Table 6). Different symbols cor- 
respond to different redshift ranges: squares z > 2.5; triangles 
2.5 > z > 1.5; circles 1.5 > z > 0.8. Absorbers for which the 
counterpart has been confirmed by spectroscopy are shown as 
filled symbols, and are labeled. The two possible DLA absorbers 
for the quasar 0841 + 129 arc joined by a dotted line. The hatched 
area shows the densely populated region from the simulation of 
Katz et al. (1996), and the dashed lines are the boundaries en- 
closing all the points in the simulation. 



corresponding in projection to the highest column densities, 
there may be a region of star formation, of diameter ~ l/i -1 
kpc. This central source would be observed as a Lyman- 
break galaxy. The gas column density decreases outwards, 
as suggested by the anticorrelation between impact param- 
eter and Nhi- The size of the region over which the col- 
umn density is > 2 x 10 20 cm~ 2 is several kpc. The region 
of star formation would be surrounded by a zone of ionised 
gas. The Ly a photons are resonantly scattered in escaping 
to the surface of the surrounding cloud of neutral gas. The 
size of the observed region of Lya emission is larger than 
the region of star formation, but smaller than the diame- 
ter of the gas cloud. The latter indicates a preferred direc- 
tion of escape, implying that the gas resides in a flattened 
structure. In reality the gas cloud is likely to be irregular 
rather than smooth, and might contain several knots of star 
formation, surrounded by HII regions of different sizes, and 
the Lya photons will escape preferentially where the column 
density of neutral gas is lowest, possibly along a complex net- 
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Figure 6. Relative positions of galaxies in four high-rcdshift 
groups discovered in Lya searches. Each group has been rotated so 
that the principal axis is horizontal in the figure, and the groups 
have been offset from one another vertically. Note the different 
scale of the lower plot. 



work of tunnels. Merging clouds will be observed as galaxies 
with irregular structure. 



5.5 Filamentary structure in the distribution of 
galaxies at z — 3 

We have previously suggested (Paper II) that the approx- 
imately linear arrangement of SI, S2, S3, as well as other 
known groups at high-redshift, may correspond to the fila- 
mentary arrangements of galaxies and galaxy sub-units that 
are found in computer simulations of the high-redshift uni- 
verse. Figure 6 summarises the observational situation, and 
shows the spatial arrangement of four high-redshift groups 
of galaxies discovered in Lya searches (from top to bot- 
tom: Francis et al., 1995, considering the Lya absorber as 
a galaxy; Pascarelle et al., 1996; Paper I; Le Fevre et al., 
1996). All groups have been rotated to a horizontal base- 
line. Clearly each of the structures is elongated. 

To quantify the degree of alignment of each group we 
have measured the smallest internal angle Sl of the trian- 
gle defined by the three objects in each group (for the PKS 
0528-250 field at z=2.81 we used the position of SI for the 
DLA absorber, rather than the position of the quasar sight- 
line through the absorber). The measured angles are, from 
top to bottom in the figure, S L = 7.0°, 1.0°, 5.8°, 2.1°. The 
average value of the angle Sl = 4.0° is certainly very sugges- 
tive of nlamentariness in the distribution of galaxies at high 
redshift. It is also notable that the sizes of the groups range 
over an order of magnitude. In the computer simulations fil- 
aments are seen at all scales (e.g. Evrard et al. 1994). We do 
not suggest that every high-redshift group of three or four 
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galaxies discovered will exhibit a similar degree of align- 
ment, but rather that filamentary networks akin to those 
seen in the simulations may become visible as larger sam- 
ples of high-redshift galaxies with denser sampling become 
available. The phenomenon might provide a useful discrim- 
inant for models of structure formation. 
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